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Abstract. We present a reanalysis of the VIRMOS-Descart weak lensing data, with a particular focus on different 

corrections for the variation of the point spread function anisotropy (PSF) across the CCDs. We show that the 

small scale systematics can be minimised, and eventually suppressed, using the B mode (curled shear component) 

VQ ■ measured in the corrected stars and galaxies. Updated cosmological constraints are obtained, free of systematics 

f^ ' caused by PSF anisotropy. To facilitate general use of our results, we provide the two-points statistics data points 

^^ , with their covariance matrices up to a scale of one degree. For the normalisation of the mass power spectrum 

f^ ■ we obtain as = (0.83 ± 0.07) {flM /0.3)~ ' . The shape parameter F was marginalised over F G [0.1, 0.3] and the 

mean source redshift Zs over [0.8, 1.0]. The latter is consistent with recent photometric redshifts obtained for the 

VIRMOS data and the preliminary spectroscopic redshifts from the VIRMOS-VVDS survey. The quoted 68% 

contour level includes all identified sources of error. We discuss the possible sources of residual contamination in 

O \ this result: the effect of the non- linear mass power spectrum and remaining issues concerning the PSF correction. 

Our result is compared with the first release of the Wilkinson Microwave Anisotropy Probe data. It is found that 

^ ' Cold Dark Matter models with a power law primordial power spectrum and high matter density Qm > 0.5 are 

excluded at 3- a. 
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1. Introduction 



200(i iKaiser et al. I I2OO0I: IVan Waerbeke at al. I I2OO0I: 
Wittman at al. Il2000(l and of tha mora numerous recent 



Weak gravitational lensmg by large scale structure (i.e., ^^^^^ ^^^^ discussions of the strength and weaknesses of 

cosmic shear) directly probes the matter distribution and ^^^j^ j^^gj^g^ h ave been recently reviewed by a number 

complements studies of the local universe based on tha pf pa pers fe.g.. IVan Waerbeke fc MellieJl200a iR.efrerieil 

hght (galaxy) distribution alone. Its sensitivity to the |2003t lHoekstrall2003a^ 

dark matter power spectrum and the geometry of the r^ ■ ^ ^ ,• ^ ■ ^^ ^ 

,. ,.;ni 1- ,-1 Cosmic shear studies are now entering the second 

universe are key ingredients tor breaking cosmological pa- , 1 ,, ^ , t^ tt •• 

, . ■ , , .,1 ,, • , phase as new surveys, such as tha Lanada-Franca-Hawaii 

rameters degeneracies associated with othe r experiments rn 1 t o /r-mTTmr om 

„ . ,,. -„ , , IT — '^^ — : — 71 lelescope Legacy burvey (CtHlLb) aim to cover areas 

such as Cosmic Microwave Backgro und Hapergel et alJ r 1 ^1 ^ i r -^ 1 1 ^i 

^^^^^^ ^^ ^^^^^'^■^^■^^'^^^^■^^^^r r\T ab"\r Thar dvo an rivHov r\T m ci rrmrnHo dTrroT" Than r'nwOTrr 



„"„„i, , ; , m T Jt-.- r 1 1 lir,^ol of sky that are an order of magnitude larger than current 

20031) and the Type la supernovae llRiess at alJ Il998t , , , rr., , ^- 



TT^ — ~ ; — , ijJnrk mi • n 1 T 7- f data sets. These second generation surveys can provide 

Perlmutter et al.ll999l) . Ihis powerful application has mo- , . ; r , • i ; •,, 
constraints on a range of cosmological parameters, with a 



tivated a first generation of surveys and the field has ..,,,. , , , , ^n ^n 

, ■ ■„ 7 ., r i-;i ;r precision that is comparable to current CMB experiments, 

made sigmhcant strides forward m the past t our years, ^t ^i i • ^i i -i-^ ^ ^ i i 

. ?, , . r ; , n ; ; Ti ' '~ , I However, the success hinges on the ability to control and 

A compilation of the nrst measurements (.Bacon et al. I ^„ , ,. , , ,. tii_li 

' correct tor observational systematics. Indeed, there are 

'Z Z £f ~Z , , u 1 **• c several technical issues remaining, concerning the point 

bend offprint requests to: waerbekeQiap.tr °' or 

* Based on observations obtained at the Canada-France- spread function (PSF) correction, the galaxy shape anal- 

Hawau Telescope (CFHT) which is operated by the National ysis and the prediction of the non-lmear power spectrum. 

Research Council of Canada (NRCC), the Institut des Sciences Whether this is possible is still debated, 
de I'Univers (INSU) of the Centre National de la Recherche 



Scientifique (CNRS) and the University of Hawaii (UH) ^ |http : //www . of ht . hawaii . edu/Science/CFHLS/| 



Van Waerbeke et al.: Results from the VIRMOS-Descart Survey 



One of the best ways to proceed is to learn from 
experience obtained from the first generation surveys. 
The VIRMOS-Descart survey that will be discussed 
below provides an excellent means to do so. It has 
already produced several re sults of c qsmological in- 
terest ([Van Waerbeke et all 12002 iHoekstra et all 
2002t Ber nardeau. Mellier fc van Waerbeke '2002'; 
Pen, van Waerbeke fc Mellier 2002: Pen et al. 2003b) . 



but showed a non-zero B-mode signal at all scales, which 
has limited its use for precise cosmological applications. 
This systematic residual signal is indeed an important 
noise contribution. It is therefore timely to understand 
the cause of B-mode generation in cosmic shear data and 
to provide better and more reliable tools to handle them. 

The purpose of this paper is to put the various is- 
sues mentioned above into perspective, and discuss their 
relevance to future surveys. We do so by presenting a re- 
analysis of the VIRMOS-Descart data, taking advantage 
of the latest improvements in PSF correction and a better 
estimate of the redshift distribution of sources based on 
recent spectroscopic data in the VIRMOS-Descart fields 
obtained by the VIRMOS-VLT Deep Survey (VVDS)2. 
Tables providing results of the best-corrected shear sig- 
nals as function of angular scale as well as errors are also 
given in this paper to help further joint analyses using 
complementary cosmological data sets. 

This paper is organised as follows: Section 2 defines 
the notation and introduces the statistical quantities to 
be used. Section 3 provides a summary of the VIRMOS 
survey. In Section 4 we discuss the impact of PSF correc- 
tion on the cosmic shear signal. Sections 5 and 6 show the 
new results from the VIRMOS-Descart survey. 

2. Theory 

For self-consistency and in order to define the notations, 
we briefly outline the revelant quantities used in this 
paper. More detailed calculations and explanations can 
be found in re c ent r eviews. Following the notation in 
ISchneider et alJ l|1998|) . we define the power spectrum of 
the convergence k as 
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where fxiw) is the comoving angular diameter distance 
out to a distance w {wh is the horizon distance), and 
n{w{z)) is the redshift distribution of the sources given in 
Ea. H12l) . P3D{k) is the 3-dimensional mass power spec- 
trum (computed from a non -linear estimation of the 
dark matter cluste ring, see e.g. lPeacock & Dodds (1996,); 
ISmith et alJ (j2003r l and k is the 2-dimensional wave vec- 
tor perpendicular to the line-of-sight. The top- hat shear 



variance (computed using a smoothing window of radius 
Oc) and the shear correlation function can be written as 
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These statistics could be derived directly from the elliptic- 
ity of the galaxies, but they do not provide a quantitative 
analysis of residual systematics. A more useful measure- 
ment of the signal and the systematics is obtained from 
the splitting the signal into its curl-free {E mode) and 
curl {B mode) components respectively. This method has 
been advocated before to help the measurement of the in- 
trinsic alignment contarnin ation in the weak lensing signal 
(|Crittenden et al.ll2001alhl ) , but it turns out to be efficient 
to measure t he residual systematics from th e PSF correc- 
tion as well (|Pen. van Waerbeke fc Mellier l Eo02'). 

The E and B modes derived from the shape of galaxies 
are unambiguously defined only for the so-called aperture 
mass variance M^ , which is a weighted shear variance 
within a cell of radius Oc- The cell itself is defined as a 
compensated filter U{9), such that a constant convergence 
K gives Map = 0: 



A4p = 



d^OKiO) U{9). 



(4) 



(Map) can be calculated directly from the shear 7 without 
the need for a mass reconstruction. It can be rewritten as 
a function of the shear if we express 7 = (7t,7r) in the 
local frame of the line connecting the aperture center to 
the gala xy. Man can therefore be expressed as function of 
7t only ()Miralda Escudelll99lt lKaiseilll992|) : 



Map- 



d207,(0) Q{e), 



(5) 



where the filter Q{6) is given from U{9): 

Q(^) = ^/' de'e'uie')-u{9) (6) 

The aperture mass variance is related to the convergence 
power spectrum (Eq^ by: 



288 r°° dfc 
^cJo ^ 



(Mlp) = ^ I ^P.{k)[Mk0c)? 



(7) 



The i?-mode is obtained by replacing jt with 7^ in 
Eq.(|SJ), and the B mode variance is denoted {M^ )±. The 
choice of U is arbitrary at this point, provided it is a 
compensated fil t er. In this paper, U{6) is chosen as in 
ISchneider et all (IT99^. 



http://www.astrsp-nirs.fr/virmos/vvds.htm 
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The aperture mass is insensitive to the mass sheet de- 
generacy, and therefore it provides an unambiguous spht- 
ting of the E and B modes. The drawback is that aper- 
ture mass is a much better estimate of the small scale 
power than the large scale power. It is obvious from the 
function J4^{k6c) in Eq.Q which peaks at k9c ^ 5, and 
essentially, all scales larger than a fifth of the largest sur- 
vey scale remain unaccessible to Af ap . The large-scale part 
of the lensing signal is lost by Map, while the remain- 
ing small-scale fraction is difhcult to interpret because the 
strongly non-linear power i s difficult to predict accurately 
l|Van Waerbeke et allEoO^ . 

It is therefore preferable to use the shear correlation 
function which is a much deeper probe of the linear regime. 
The E and B modes can also be measured separately 
from the shear variance and the shear correlation functions 
(Eq[21l2Il- However, in contrast with the Map statistics, the 
separation of the tw o modes can only be don e up to an 
integration constant (jCrittenden et al.ll2001bj) . which de- 
pends on the extrapolated signal outside the measurement 
range, either at small or large scale. 

An alternative, which does not require the knowledge 
of the signal outside the measured range, is to use the 
aperture mass B mode to calibrate the shear correlation 
function B mode. An range A9 of angular scales where 
(Map {A9))± ^ would guarantee that the B mode of the 
shear correlation function should be zero as well (within 
the error bars), at angular scales ^ 5A9. The practical 
implementation of this calibration scheme is discussed in 
Section 5. 

The E and B modes of the shear top-hat variance and 
correlation function are accessible from the following shear 
correlation function ^+ and ^_ : 



Mr) = ^M0ht{0 + r)) + {jriehri0 + r)). 

?_(r) = M9ht{9 + r)) - (7r(e)7.(^ + 0), (9) 

where 7t and 7^ are the tangential and radial projec- 
tions of the shear on the local frame joining two ga laxies 
separat ed from a distance r. Following .Crittenden et alJ 
I 2nniah . we define 
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The E and B shear correlation functions are given by 
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In this paper, we use Ea. pi|) as the cosmic shear signal for 
constraining the cosmological model. The aperture mass 
(Eq[7I) is used only to calibrate the B mode of the shear 
correlation function, at scales where {M^ )± is measured 
to be zero. 



3. The VIRMOS data set, redshift distribution 
and statistical estimators 

We use the observations carried out within the VIRMOS- 
DESCART project 3 by the VIRMOS ^ i maging and 
spectroscopic survey (|Le Fevre et al.l l2004|) . The data 
cover an effective area of 8.5 sq.deg. (12 sq.deg. be- 
fore masking) in the I-band, with a limiting magnitude 
Iar = 24.5. Technical deta ils of the data set are given i n 



IVan Waerbeke et alJ l)200l[) and lMcCracken et all | 20& 



We applied a bright magnitude cut at Iab = 21 in or- 
der to better exclude foreground objects from the source 
galaxies. 

As we already described in I Van Waerbeke et alJ 
1)20021) . we use an estimate of the source redshift distri- 
bution based on photometric redshifts measured in joint 
VLT and HDF data. The source redshift distribution is 
parametrised as 



n{z) 



P 



(T) 



-, exp 



(12) 



with a = 2, /3 = 1.2 and Zg = 0.44. For these values of a 
and j3, the mean redshift is Zg w 2.1 Zg and the median 
redshift is w 1.9 Zg. These values are given by the first 
VVDS data sets, and provide an accurate calibration of 
the source redshift distribution used in VIRMOS-Descart. 
It agrees nicely with spectroscopic redshift measurements 
in the VVDS survey (a preliminary release can be found 
in iLe Fevre et al.l 12003). The uncertainty in the source 
redshift distribution is discussed in Section 5. 

The galaxy shapes are meas ured and ana l ysed w ith 
IMCAT, which is described in iKaiser et all l|l995l) to 
which we refer for technical details. This technique allows 
us to measure, for each galaxy i, an ellipticity e{Oi) = 
(61,62) and a weight Wi for each galaxy and star in the 
data. The ellipticity is given by a weighted second order 
moment of the light distribution f{0) of the galaxy which 
is an unbiased estimate of the shear 'y{9i). The signals 
measured from the galaxy shapes are the binned tangen- 
tial and radial shear correlation functions. These are given 
by a sum over galaxy pairs {6i, 9j) 



Ur) 



^rr{r) = 



'^WiWjet{Oi) ■ et{Oj) 

'^WiWjeriOi) ■ er{9j) 
hi 

i,3 



(13) 



where r = \9i — 9j\, and (6t,e,.) are the tangential and 
radial ellipticities defined in the frame of the line connect- 
ing a pair of galaxies. The weights Wi are computed from 

^ http://terapix.iap.fr/DESCART 
^ http://www.astrsp-mrs.fr 
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the intrinsic ellipticity variance crl and the r.m.s. of the 
elUpt icity PSF correction a^ (f or the details, see Section 
3 in IVan Waerbeke et allEoog) . ^ 



4. Updated analysis of the VIRMOS-Descart data 

4.1. Previous sources of B-mode 

IVan Waerbeke et alJ l)2002|) discussed the presence of 
residual systematics, and attributed the B mode to an 
imperfect correction for PSF anisotropy The presence of 
this B mode limited the precision with which the lensing 
signal could be determined. The -B-mode was comparable 
to the la error. As it was not clear how to reduce it and 
since there was also no obvi ous way to quantify the con- 
tamination on the E mode, IVan Waerbeke et al.l l|2002() 
did not consider any further improvement. Rather, they 
included the systematics in the final error budget, lead- 
ing to a slightly biased normalisation of the mass power 
sp ectrum towards h igh values. 

iHoekstral ()2004 studied the effects of PSF anisotropy 
using CFHT observations of fields with a high number 
density of stars. This study revealed that the second or- 
der polynomial model for the PSF anisotropy used in 
IVan Waerbeke et al.l l|2002|) does not ac curately describe 
the spatial variation of the systematics. iHoekst ra ( 2004) 
showed that there are rapid changes in the PSF anisotropy 
at the edges of the CFH12k mosaic. In the actual VIRMOS 
images, the small scale PSF correction is made difficult by 
the lack of stars. The mean star separation at high galactic 
latitude is slightly less than 2'. Therefore the choice of the 
correct PSF correction model must play an important r ole 
for the small scale cosmic shear signal l)Hoekstrall2004l) . 

In order to estimate the effect of model choices, we 
tested various PSF models for the spatial variation of the 
PSF anisotropy. The ellipticities of the observed stars are 
modeled using generic functions: 



2x10-" 
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1 + E;,m ^'™ fim{x,y)' 



(14) 



where the fim{x,y) are the generic formal polynomial 
function terms listed in Table ^ and the aim and bim are 
their coefficients, as derived from a least squares fit of 
the model to the sample of stars selected from the data. 
Tabic n lists six of the PSF models that were considered. 
Models 1 — 5 are simple polynomials (i.e., bim = 0), with 
higher orders in y because the CFH12k chips are longer 
in that direction and because the rapid change in PSF 
anisotropy is a strong function of y. Model 6 is a "tai- 
lor" made mode l, based on the study of dense star fields 
llHoekstralbOO^ . 

For each model, we measured the ellipticity correlation 
functions from the stars (after correction) and computed 



In practice, r is the averaged angular distance computed 
inside very tiny bin size (« 3 pixels only). Note that the r' 
value of the shear statistics is different and correspond to the 
upper limit: r < r' 




-1x10- 





IxlO-n 



10-" 



■1x10- 




*V^ 



10 



10 



6 rarcmin] 

ap L J 

Fig. 1. The E (black area) and _B- modes (light/green 
area) aperture mass statistics computed from the residual 
shapes of the stars after correcting using on of the six 
PSF models discussed in the text. The areas indicate the 
la regions around the mean. 



the corresponding E and B mode aperture mass statis- 
tics. The results are presented in Figure Q] We find sig- 
nificant deviations from E = and B — at small scale 
(< 10') for all models. Note that these results do not tell 
us about the amount of residual systematics in the galaxy 
sign al, although th ey can be used to make an estimate 
(see lHoekstrall2004J) . 

In the absence of more detailed information. Figure Q 
provides a very useful way to choose which of the PSF cor- 
rection models works best. The rational function (model 6) 
shows the smallest contamination on all scales. In partic- 
ular the deviation from zero starts at a smaller scale than 
the other models. FigureHdemonstrates there is still room 
left for an even better correction at small scales, but the 
lack of suitable stars makes this a very challenging issue. 

These new insights in the correction for PSF 
anisotropy have led to a significant reduction in the level 
of systematics of the VIRMOS-Descart cosmic shear sig- 
nal, as we will demonstrate below. Several other s ources 
of the B mode observed by IVan Waerbeke et al.l 1)200211 
were identified as well. The systematics that have been 
corrected for in the revised analysis presented here are: 

1. The small scale B mode was generated by an imper- 
fect PSF correction mode l as described above. As sug- 
gested in IHoekstral l|2004l) a proper correction for PSF 
anisotropy lowered the E mode signal on these scales. 
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Table 1. Different Point Spread Function parametric 
models used to correct the star anisotropy in the VIRMOS 
data. The table lists the basis of functions of the mod- 
els. Model 6 is a rational function, the other models have 
bim = 0. 



model # 


Correction Model 


1 
2 
3 

4 
5 
6 


{l,x,y,x'^,xy,y'^,x^,y^,y*') 

{l,x,y,x^,xy,y^,y^) 

(1, X, y, x^,xy, y^,y^,y'^, y^,x-^,x^) 

{l,x,y,x^,xy,y'',x\y^) 

(1, X, y, x'^,xy, y^,y'^, y )/{x, y) 
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Fig. 2. top panel: E (filled points) and B (open points) 
mode aperture mass measurements as measured from the 
previous analysis. It uses model-l (see Table ^ to de- 
scribe the PSF anisotropy with the wrong centroid esti- 
mate and no sigma clipping on the shape parameters of 
the selected stars, bottom panel: results using model 6 and 
the improved catalogues. Because of a different object se- 
lection in the new version, we ensured that we only used 
objects that are in common to produce this Figure. The 
new catalogue contains more galaxies. 



The large scale constant B mode was caused by an im- 
perfect calculation of the centroid of the objects. This 
problem was the most difficult to find because it was 
caused by a subtle inconsistency between the centroid 
estimate (based on SExtractor) and the shape calcu- 
lation (based on IMC AT). On one of the VIRMOS 
fields (FIO), SExtractor systematically computed a bi- 
ased centroid which was the cause of the large scale 
bias in shape measurement generating a B mode. This 



problem was cured by measuring the centroids using 
aperture filters instead of isophotal apertures. 
3. The third important source of systematics arose from 
spurious ellipticies for some of the selected stars. Some 
of the stars survived the selection criteria based on a 
sigma clipping applied to their corrected ellipticity. It 
turned out that some stars passed this selection cri- 
teria, but their smear polarisability tensor was signifi- 
cantly d ifferent from the bu lk of the stars (for a defini- 
tion see iKaiser et al.ll995l) . This could happen when a 
star has a small very close companion which would con- 
taminate the derivative of the profile much more than 
its ellipticity. The application of a 3 — (t clipping to all 
of the star parameters, corrected ellipticity, smear po- 
larisability and the stellar part of the preseeing shear 
polarisability remove all the remaining unwanted stars 
which contributed to the residual B mode. 

Figure |21 demonstrates how the results have changed 
from the early VIRMOS analysis (top panel) and the cur- 
rent analysis (bottom). The top panel shows clearly the 
small and large scale residual B mode as measured before 
and a "bump" in the E mode at '--^ 5 arcminutes. The bot- 
tom panel shows the improvement after the main residual 
systematics have been corrected for. It also shows that 
the aperture mass is an effective statistic to identify and 
quantify residual systematics. 



4.2. Residual systematics 

The results of the new analysis of the VIRMOS-Descart 
are presented in Figure |3| To derive these measurements, 
all sources of systematics discussed in the previous sec- 
tion have been accounted for, and model 6 was used to 
characterize the spatial variation of the PSF anisotropy. 
Figure Ol shows the resulting ellipticity correlation func- 
tion as a function of angular scale. The measurements are 
tabulated in TableEl As explained in Section 2, in order to 
separate the signal into an E and B mode for the top-hat 
shear variance and shear correlation function, we have to 
define a zero-point for the B mode. Based on the results 
for the aperture mass statistic (panel b) the correlation 
functions are scaled such that the mean B mode is zero 
on scales 2 — 10' according to the vanishing B mode for 
the aperture mass at 6* = 10 — 50'. The error bars reflect 
only the statistical error bars. There is no cosmic variance 
on the B mode, and consequently the statistical error is 
an unbiased estimate of the B mode noise. 

Panels (b) and (c) show the aperture mass (M^^) and 
the shear top-hat variance (7^), respectively. The mea- 
surements of both statistics are listed in Table |21 The 
small scale (Mf ) B mode is consistent with no signal 
(i.e., no systematics). However, one has to bear in mind 
that different PSF correction models produce fluctuations 
up to 10% in the amplitude of the aperture mass below 
1'. On the other hand, the aperture mass signal is robust 
against changes in the adopted PSF correction model for 
the range of scales 10 — 50' which corresponds to a phys- 
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Fig. 3. Panel (a): E and B mode shear correlation function (filled and open points, respectively) as measured from 
the VIRMOS-Dcscart data; panel (h): E and B mode measurements for the aperture mass statistic; panel (c): as 
before but now for the top-hat variance. For all three statistics we observe the signal is consistent with no i?-mode on 
large scales, and a hint of minor contamination on small scales. 



ical angular scale range of 2 — 10'. This means that the 
smallest scales in the aperture mass statistics give robust 
cosmological signal only if the PSF variation is properly 
removed. 

The amount of residual systematics left in the signal 
due to imperfect PSF correction can be estimated from 
the correlation between the uncorrected stars an d the cor- 
rected galaxies. Such an estimator was defined in lHevmanj 
(|2f 



CsYS 



{e*l? 



(e'e" 



(15) 



where e* is the star ellipticity before PSF correction, and 7 
is the shear estimate of the galaxies, ^gys is renormalised 
by the star ellipticity correlation function (e*e*), which 
makes it directly comparable to the signal (7(^)7(6' + r)), 
and not dependent on |e*|. 

Figure il shows ^sys for the VIRMOS data. The resid- 
ual PSF contamination is consistent with no signal except 
maybe for one point around 6 ~ 22arcmin, which is at 
2.4(7. It demonstrates that the anisotropy of the PSF has 
been almost completely removed. However, this plot pro- 
vides no information about the accuracy of the isotropic 
PSF correc tion (the "pre-see i ng" sh ear polarisability). As 
outhned in iHirata fc SeliakI l|2003|) . the accuracy of the 
isotropic correction of the PSF is still somewhat uncertain, 
essentially because there is yet no direct way to obtain 
a perfect calibration of the shear amplitude to compen- 
sate for the PSF smear ing. One s hould emphasize that 
lErben et al. (2 001)) and iBacon et a l. ( 2GGI5) used simu- 
lated images to demonstrate that the shape mea surement 
method used in this paper iJKaiser et al.Ml995|) is accu- 



Table 2. E and B modes of the shear correlation func- 
tion. The error is statistical only (see Section 6 for the 
covariance error). 



e 


iE 


Cb 


5^ 


0.780000 


0.000200166 


5.434006-06 


3.330006-05 


1.12000 


0.000122166 


5.534006-06 


3.250006-05 


1.61000 


0.000109166 


1.743406-05 


1.850006-05 


2.45000 


9.75660e-05 


2.183406-05 


1.320006-05 


3.42000 


4.56660e-05 


-4.60600e-06 


1.140006-05 


4.39000 


4.43660e-05 


9.034006-06 


1.030006-05 


5.37000 


2.80660e-05 


8.140006-07 


9.460006-06 


6.83000 


3.40660e-05 


4.734006-06 


6.290006-06 


8.79000 


2.48660e-05 


2.034006-06 


5.710006-06 


12.2000 


2.79660e-05 


-2.24600e-06 


3.380006-06 


17.0800 


1.39860e-05 


7.534006-06 


2.980006-06 


21.9600 


1.46060e-05 


3.040006-07 


2.760006-06 


27.0100 


1.18960e-05 


1.934006-06 


2.510006-06 


32.3300 


1.10960e-05 


2.034006-06 


2.320006-06 


37.5400 


7.80100e-06 


-2.06000e-07 


2.310006-06 


42.5400 


8.22120e-06 


-7.46000e-07 


2.230006-06 



rate to 5% or better at recovering the correct shear am- 
plitude. The worst calibration was obtained for the most 
ellipticial galaxies. Hence, we expect an accuracy due to 
isotropic calibration smaller than 5% in this study, which 
is still smaller than the statistical error. The PSF anal- 
ysis technique based on reducing the B mode amplitude 
shows how it is possible to control the PSF anisotropy 
correction, however the isotropic correction still has to be 
checked using simulated images. 
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Table 3. Values of the aperture mass and shear top-hat variance as a function of scale. The B mode is rescaled to 
zero for the later, and the error is statistical only. 



e 


(M^P> 


(A4^p)x 


<5(A/i) 


(7l> 


(Ib) 


S{l') 


0.59 


0.1956E-04 


-.2560E-04 


0.1571E-04 


0.000246658 


5.66021e-05 


2.187006-05 


0.98 


0.1806E-04 


0.8109E-05 


0.9666E-05 


0.000183258 


3.11921e-05 


1.419006-05 


1.27 


0.2318E-04 


0.9647E-05 


0.7599E-05 


0.000158658 


2.56821e-05 


1.147006-05 


1.97 


0.1703E-04 


0.2864E-05 


0.5083E-05 


0.000122858 


1.90721e-05 


8.005006-06 


2.94 


0.1378E-04 


0.2319E-05 


0.3512E-05 


9.088796-05 


1.26021e-05 


5.699006-06 


3.92 


0.1395E-04 


0.2413E-05 


0.2704E-05 


7.15379e-05 


8.56213e-06 


4.457006-06 


4.89 


0.1357E-04 


0.2952E-05 


0.2217E-05 


5.922796-05 


6.85213e-06 


3.688006-06 


5.87 


0.1250E-04 


0.2675E-05 


0.1884E-05 


5.090796-05 


5.65212e-06 


3.155006-06 


7.82 


O.llOOE-04 


0.1198E-05 


0.1465E-05 


4.186796-05 


3.34212e-06 


2.472006-06 


9.76 


0.8562E-05 


0.1469E-05 


0.1210E-05 


3.633796-05 


2.52213e-06 


2.049006-06 


14.65 


0.4570E-05 


0.1084E-05 


0.8642E-06 


2.724796-05 


2.59212e-06 


1.467006-06 


19.52 


0.3899E-05 


0.1008E-06 


0.6924E-06 


2.185796-05 


2.45213e-06 


1.171006-06 


24.41 


0.3594E-05 


-.2831E-06 


0.5858E-06 


1.835696-05 


1.92213e-06 


9.881006-07 


29.62 


0.3503E-05 


-.4899E-06 


0.5065E-06 


1.615396-05 


7.78125e-07 


8.553006-07 


35.05 


0.3537E-05 


-.3382E-06 


0.4477E-06 


1.503296-05 


-3.318766-07 


7.568006-07 


40.05 


0.3197E-05 


0.4226E-07 


0.4088E-06 


1.435696-05 


-1.409886-06 


6.896006-07 


45.05 


0.2635E-05 


0.2219E-06 


0.3804E-06 


1.375096-05 


-2.507886-06 


6.382006-07 


49.95 


0.1956E-05 


0.3912E-06 


0.3599E-06 


1.301896-05 


-3.494886-06 


5.992006-07 



10- 



10-" - 



?-- 



10- 




10-8 

10 
0[arcmin] 

Fig. 4. The shades (red) area corresponds to the la area 
around the mean shear correlation function {'y{r)^{9 -\- 
r)) measured from the VIRMOS-Descart data. The points 
with error bars show the residual systematics correlation 
function ^SYS as defined in the text. The indicated error 
bars are those of the shear correlation function, and are 
displayed to show the level of systematics with respect to 
the 1 — (7 statistical error. 



5. Parameter estimation 

5.1. Likelihood calculation 



Following IVan Waerbeke et al.l l)2002|) . we investigate a 
four dimensional parameter space defined by the mean 
matter density Qm, the normalisation of the power spec- 
trum cTg, the shape parameter T and the redshift of the 
sources parameterized by Zs (see Eg. lfT^ 'l. We adopt a flat 
cosmology, unless stated otherwise. The default priors are 
taken to be fiM e [0.1,1], erg G [0.5, 1.3], T e [0.1,0.5] and 
Zs £ [0.38, 0.48]. The latter corresponds to a mean redshift 
betwe en 0.8 and 1.0, in agreement with iLe Fevre et alJ 

(Hooi). 

The data vector ^i is the shear correlation function 
£,E as a function of scale, as listed in Table El This table 
also shows that the residual systematic £,b is negligible. If 
we take mi{flM, as,T, Zg) to be the model prediction, the 
likelihood function of the data is given by 

^ ^ (2^)n|c|l/^ '"P ^^^' - "^^)*^"'(^' - ^^^^^ ' (16) 

where n = 16 is the number of angular scale bins and C 
is the 16 X 16 covariance matrix. 



Cij = (te ~m.if{^j -nij)). 



(17) 



C can be decomposed as C = C„ + C^, where C„ is the 
statistical noise and Cs the cosmic variance covariance 
matrix. C„ is diagonal and given in Table El(last columns). 
The m atrix C., is computed according to lSchneider et alJ 
l|2002|) . assuming an effective survey area of 8.5 square de- 
grees, a number density of galaxies Ugai = 15 per arcmin^, 
and an intrinsic elliptic ity dispersion a^ = 0.44. These 



numbers can be found in lVan Waerbeke et al.l 1 200211 . The 
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cosmic variance is comp uted assuming a Gaussia n statis- 
tic, but as shown in iVan Waerbeke et al.1 l)2002J) . this is 
not an issue in our case: the non gaussian contribution ap- 
pears only at small scale, where the noise is still dominated 
by the statistical error. We need a fiducial model to com- 
pute the matrix components, which i s taken as the best fi t 
model following the WMAP results ijSoergel et alJl2003|) . 
It corresponds to the choice Hm = 0.3, A = 0.7, erg = 0.85, 
r = 0.21 (the reduced Hubble constant is h=0.7). We take 
a mean source redshift for the fiducial model of 0.9, which 
corresponds to Zg — 0.43. Table0]gives the full covariance 
matrix C used in this paper. The pure cosmic variance 
part of the matrix (Cj) can easily be extracted by sub- 
tracting the statistical noise contribution given in Tablc|21 
The B mode is calibrated by marginalizing around B = 
within the la interval. 

The results of the maximum likelihood analysis are 
presented in Figure [Sj which shows the constraints in 
the Qm-cs pl an for the non-li near prediction schemes 
described in iSmith et alJ 1120031). We not e that we find 
similar results for the IPeacock fc Doddsl l|l99(il) model, 
although the best fit cosmologies differ somewhat (e.g., 
see Figure El). We marginalised over F e [0.1,0.3] and 
Zg e [0.4,0.48], assumi ng a fiat probability distribution in 
these intervals. For the lSmith et al.l 1 20031) model we find 



as = (0.83 ± 0.07) ( ^ 



-0.49 



(18) 



For the IPeacock fc Doddsl l|1996|) model we find a some- 
what higher normalisation, given by 



as = (0.86 ± 0.08) 



"(13 



-0.49 



(19) 



The difference between these two non-linear models 
corresponds to a 4% discrepancy, and it shows the lim- 
itation in the cosmic shear signal interpre tation due to 
non-linear modeling. As the ISmith et al.l l|2003ll is the 
most recent, and more extensive analysis of the non- 
linear power spectrum, we prefer the results of this model. 
The updated results are also in excellent agreement with 
results obtained from the Red-sequence Cluster Survey 
(Hoekstra ct al. 200i), who found as = 0.86t!]:g^ for 
^m = 0.3. They are also in agreement wit h the revisited 
CTIO lensing survey ((.Tarvis fc Jain II2004I) . 

Figure El shows the probability distribution functions 
for as for the two models, adopting flm = 0.3. Despite the 
different maximum likelihood shifts between the models, 
the curves are well fitted by a gaussian distribution as 
shown by the dashed line. 

It is interesting to note that the new PSF anisotropy 
co rrection drops _ the res ulting as only by 10% as compared 
to lVan Waerbeke et alJl)200l|) . much less than what would 
be expected from the change in amplitude in the lensing 
signal. However, in I Van Waerbeke et al. (200l|). the source 
redshift distribution was given by the model described in 
IWilson et al. I l|200lh (their figure 4), which turns out to 




Fig. 5. Constraints on as and flM for a flat CDM 
Universe, adoptin g the non-linear prediction given by 
ISmith et al.1 Ij2003r . The contours correspond to the 0.68, 
0.95 and 0.999 confidence levels. The solid dark vertical 
line indicates Q = 0.3 and the horizontal line indicates the 
maximum of the likelihood for as sd, fl — 0.3. The models 
are marginalised over F G [0.1,0.3] and z^ e [0.4,0.48]. 



give too much weight to high redshift galaxies compared 
to w hat we know a bout the VIRMOS galaxy sample to- 
day IILe Fevre et al. 2003. 2004). . Hence the residual sys- 
tematics in IVan Waerbeke et al.l (|200ll) were partly com- 
pensated by this biased redshift distribution model. For 
the present analysis we use more recent results from a 
comparison with photomet ric and spectroscopic redshifts 
(|Le Fevre et al.l2003ll20ol . This changes the redshift dis- 
tribution and significantly reduces the associated uncer- 
tainties. 



5.2. Combination with \NMAP 

Figure m indicates that a high density Universe SI a/ = 1 is 
excluded at the 2-cr level. This result, however, depends on 
the prior used for the shape parameter F r^ ^uh (which is 
F S [0.1, 0.3]), and therefore on the Hubble constant (and 
the baryon abundance to a less extent). Using weak lens- 
ing data alone, it is still possible to accommodate ^u — 1, 
provided F is at least as large as 0.5, and we cannot 
exclude the high density {Hm > 0.9) alternative to th e 
concordance model, proposed bv lBlanchard et al.1 l|2003l) . 
The latter model only works if the Hubble constant is low 
(~ 0.5 , which is not in agreement with IFreedman c t alJ 
Ij200ll) 'l and if the Universe expansion rate is not accel- 
erat ing as sug gested by t he type la sup ernovae results 
(Rie ss et al.l ll998: Pcrlmutter et al.l ll999l) . 

From the cosmic shear point of view, such a model 
would require a low normalisation as — 0.4, because for 
higher values of F, the banana shape contour in the ^m — 
as plane (Figure EJ is tilted clockwise. For example, the 
combination of F = 0.5 and as = 0.5 is excluded at more 
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Table 4. Full covariance matrix C, as discussed in the text, Section 6. Units are 10^^". The scales 9i correspond to 
those given in Table |2| 
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6*5 
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^12 
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01 


14.48 


2.15 


1.85 


1.55 
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1.20 
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0.92 


0.80 


0.70 


0.62 


0.56 
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02 
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0.49 
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6i 
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0.49 


0.44 
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Fig. 6. a^ probability distribution functions for rtf,[ = 
0.3. Dashed curve is for |,Pcacock & Dodds (1996) (high 
normalisation), and the solid curve for lSmith et alJ <|2003(1 
(low normalisation). The superimposed dotted lines indi- 
cate a Gaussian fit to the likelihood function. 



than 3-a. This is inconsistent with WMAP IjSpergel et alJ 
|200,1J) which gives us = 0.9 ±0.1 for a flat and adiabatic 
Cold Dark Matter model assuming a power law primordial 
power spectrum. 

Combined with our cosmic shear results, this leads to 
a rejection of the very low and high density universes, 
with a preferred value around Q.m — 0.3, as indicated by 
Figure[7| This result also demonstrates that CMB and cos- 
mic shear experiments are complementary, enabling one to 
break the Q. m — ts degeneracy without ex ternal data sets 
ijContaldi et all 120031: iTerenoet al.ll2004 . Although one 
can modify the matter content of the Universe to try to 



Fig. 7. VIRMOS-Descart data constraints in the A — J^m 
plane. The models are normalised to the WMAP results 
((78 — 0.9 ± 0.1), and we assumed the priors P = [0.1, 0.5] 
and Zs = [0.4,0.48]. 



account for a low erg at high Qm IjBlanchard et al.ll2003|) . 
the next generation of weak lensing survey will certainly 
give an unambiguous constraint on cg at low redshift re- 
gardless the value of fl]\i . 

6. Conclusion 

We have presented a revised analysis of the VIRMOS- 
Descart cosmic shear data, focussing on how to deal with 
residual systematics. The main result of this paper is a 
10% smaller mass power spectrum normalisation com- 
pared to our previous estimate. We find erg = 0.83 ± 0.07 
(for the preferred non-linear power spectrum as described 
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by Smit h et alJ 1)2003 !)'). which is at t he 1-g bo ttom end 
of the measurement given in Van Waerbeke et al. (2002) , 
where the systematics discussed in this study were stih 
present. Our better understanding of systematics allows 
us to use weaker cosmic shear signal at higher angular 
scale and therefore to be less sensitive to unknown prop- 
erties of the non-linear power spectrum. 

In this work, various sources of systematics, which were 
previously ignored, were identified and the amplitude of 
the residual B mode from the aperture mass was used to 
significantly improve the cosmological lensing signal. The 
choice of a better model to describe the spatial variation of 
PSF anisotropy at small scale proved to be critical for an 
unbiased measurement. We demonstrated that small dif- 
ferences between correction schemes can have a significant 
impact on the small scale signal. Fortunately, the residual 
systematics in the aperture mass of the stars provide an 
unambiguous method to test for residual systematics by 
selecting the least biased model. 

On large scales, a few outliers in the star catalogue 
used to correct for the effects of the PSF can introduce 
a bias even if their ellipticities appear reasonable. It is 
therefore important to reject spurious stars not only from 
their corrected ellipticity, but using extended criteria such 
as the smear polarisability and the preseeing shear polar- 
isability We also found that a large angular scale bias can 
be introduced from the use of galaxy centroids based on 
isophotal cuts. 

There is a fundamental limit to the accuracy of the 
correction at small scales, which is determined by the sam- 
pling of the PSF variation across the CCDs. Further high 
precision cosmic shear studies will require the o bservation 
of den se star fields to improve the correction (see iHoekstral 
12004 1). Fortunately, the large scale signal (> 10') appears 
unaffected by the choice of the PSF anisotropy model. 

The use of the shear correlation function is much less 
sensitive to systematics than the shear top-hat variance or 
even the aperture mass variance. This is most critical for 
the aperture mass which is sensitive to very small scales. 
On the other hand, the aperture mass B mode is a crucial 
test of the residual systematics, which has to be used to 
calibrate the PSF correction and the B mode of the shear 
correlation function. 

The prediction of the non-linear power spectrum 
results in values for as that differ by a few per- 
cent for different models. The model presented by 
Igmit li et al. ( 200^ always gi ves a smaller normalisation 
than lPeacock fc Doddsl l|199fil) . Currently, the former pro- 
vides the preferred model. We also note that these differ- 
ences are not yet an important issue for surveys of the 
size we considered here. However, it will be a challenging 
task for the forthcoming lensing surveys to improve the 
non-linear predictions to the required level of precision. 

It is still possible to accommodate a large value for 
flm from cosmic shear measurements alone. However, the 
combination of the cosmic shear signal with measurements 
of the Cosmic Microwave Background provides strong ev- 
idence for a low density universe without adding exter- 



nal data sets. To derive these conclusions we do have to 
assume a Cold Dark Matter scenario with a primordial 
power law power spectrum. Such assumptions can be re- 
laxed with the next generations of lensing surveys. 
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